Exact analytical solution of average path length for Apollonian networks 
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With the help of recursion relations derived from the self-similar structure, we obtain the exact 
solution of average path length, d t , for Apollonian networks. In contrast to the well-known numerical 
result dt oc (lniVi) 3/4 [Phys. Rev. Lett. 94, 018702 (2005)], our rigorous solution shows that the 
average path length grows logarithmically as dt oc In JVj in the infinite limit of network size N t . The 
extensive numerical calculations completely agree with our closed-form solution. 
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One of the most important properties of complex net- 
works is average path length (APL), which is the mean 
length of the shortest paths between all pairs of vertices 
(nodes) Q. Most real networks have been shown to be 
small- world or ultra small- world networks 0, S 0, 
that is, their APL d behaves a logarithmic or double log- 
arithmic scaling with the network size N: d ~ In AT [6[ 
or d ~ In In AT jjjj]. It has been established that APL 
is relevant in many fields regarding real-life networks. 
In the design or interpretation of routes in architectural 
design, signal integrity in communication networks, the 
propagation of diseases or beliefs in social networks or of 
technology in industrial networks, APL is a natural net- 
work statistic to compute and interpret. It is strongly 
believed that many processes such as routing, search- 
ing, and spreading become more efficient when APL is 
smaller. So far, much attention has been paid to the 
question of APL @, |, M, [ll|, QJ, Q]| . 

Recently, on the basis of the well-known Apollonian 
packin g [l4l|. Andrade et al. introduced Apollonian net- 
works [15] which were also proposed by Doye and Massen 
in Ref. [161 ] simultaneously. Apollonian networks belong 
to a deterministic growing type of networks, which have 
drawn much attention from the scientific communities 
and have turned out to be a useful to ol [l7l . [H, [Til [2(| 

m mmmmnEiinini, ^ ^^ Many topo ~ 

logical properties of Apollonian networks such as degree 
distribution, clustering coefficient, and correlations have 
been determined analytically [l5|, [l6| , and the effects of 
the Apollonian networks on several dynamical models 
have been intensively studied, including Ising model and 
a magnetic model [la, HI, HHi EH ■ Despite the importance 
and usefulness of the quantity APL, there is no analytical 
calculations for the APL of Apollonian networks. 

In this report, we derive an exact formula for the aver- 
age path length characterizing the Apollonian networks. 
The analytic method is based on the recursive construc- 
tion and self-similar structure of Apollonian networks. 




FIG. 1: (a) Second construction method of Apollonian net- 
work. At+i, is obtained by joining three copies of At denoted 
as A^ (tp — 1,2,3), which are connected to one another at 
the edge nodes (i.e., X, Y, Z, and O). Note that the central 
node is represented as O (marker not visible), (b) Illustra- 
tion of the recursive definition of node classification. From 



the node classification of A 



(1) A ( 2 ) 



and Ai 3 \ we can derive 
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recursively the classification of nodes in network A t +i. 



Our rigorous result shows that APL grows logarithmi- 
cally with the number of nodes. The obtained analytical 
solution modifies the previous numerical result in [To] ], 
where the authors claimed that the APL of Apollonian 
networks scales sub-logarithmically with network size. 
Our analytical technique could provide a paradigm for 
computing the APL of deterministic networks. 

The Apollonian network, denoted as At (t > 0) after 
t generations, is constructed as follows [la ]: For t = 0, 
At is a triangle. For t > 1, A t is obtained from A t ~\. 
For each of existing triangles created at step t — 1, we 
add a new vertex and join it to all the vertices of this 
triangle. Alternatively, Apollonian network can be also 
created in another method [35|, [36j . Given the generation 
t, A t +i may be obtained by joining at four edge nodes 
three copies of A t , see Fig.QJa). According to the latter 
construction algorithm, we can easily compute the total 
number of vertices of A t is Nt = (3* + 5)/2. 

Apollonian network presents the typical characteris- 
tics of real- life networks in nature and society @, l37l ]. 
It has a power-law degree distribution with exponent 
7 = 1 + In 3/ In 2 pl| [l6j], which belongs to the inter- 



2 



val between 2 and 3. For any individual vertex with 
degree k, its clustering coefficient C{k) is also approxi- 
mately inversely proportional to its degree k as C(k) = 
(6/fc) — [2/(fc — 1)]. The mean value C of clustering co- 
efficients of all vertices is very large, which asymptot- 
ically reaches a constant value 0.8284. Moreover, the 
network is small world, and its diameter, defined as the 
longest shortest path length between all pairs of vertices, 
increases logarithmically with the number of vertices [38| . 
In fact, Apollonian network can be expanded to general 
cases fl6l KM] . l38j associated with other self-similar pack- 
ings [39|]. Analogously, we can construct the stochastic 
version of these networks [HI, 0, |42j • The main topo- 
logical properties of the Apollonian type networks are 
controlled by the dimension of packings. 

After introducing the Apollonian networks, we now in- 
vestigate analytically the average path length. We rep- 
resent all the shortest path lengths of network A t as a 
matrix in which the entry g uv is the shortest path length 
from node u to v. A measure of the typical separation 
between two nodes in A t is given by the average path 
length d t defined as the mean of shortest path lengths 



over all couples of nodes, i.e., d, 



where 



N t (N t -l)/2> 

Dt = J2u veA ujtv 9 UV denotes the sum of the shortest 
path length between two nodes over all pairs. 

By the second construction, it is obvious that Apollo- 
nian network has a self-similar structure. Thus, the total 
distance D t +\ satisfies the recursion relation 
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t+i 



3A + e t 



(i) 



where 0t is the sum over all shortest path length whose 
endpoints are not in the same A t branch. The solution 
of Eq. (J) is 



D t = 3 t - 1 D 1 + (3* 



(2) 



Thus, all that is left to obtain D t is to compute m . 

The paths that contribute to t must all go through 
at least one of the four edge nodes (X, Y, Z, and O) 

at which the different At branches (A\ , A^\ A^) are 
connected. The analytical expression for t , named the 
crossing path length, can be derived as below. 

Denote 0"' /3 as the sum of all shortest paths with end- 



't 



points in Af' and A[^ . Note that 0"' p rules out the 



a,/3 



paths with endpoint at the edge of A\ a ' and Af' . For 



1 0?) 



example, each path contributed to 0j' 2 should not end 
at node O or X. Then the total sum Q t is given by Ot = 



1.2 







that 6 t 



1,3 



2,3 



By symmetry, 9 t 1,2 = 6 t 1,3 = <d z t 



2,3 
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3 0j' 2 , where &]' 2 is given by the sum 
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(3) 
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,v^X\Ju,v^O 



To calculate the crossing path length <d t ' , we classify 
nodes in network A t+ \ into seven different parts accord- 
ing to their shortest path lengths to each of the three 



vertices (i.e. X, Y, Z) of the peripheral triangle AXYZ. 
Notice that vertices X,Y, and Z themselves are not clas- 
sified into any of the seven parts represented as Pi, Pi, 
P3, P4, P5, Pq, and Py, respectively. The classification of 
nodes is shown in Fig.[ljb). For any interior node v, we 
denote the shortest path lengths from v to X, Y, Z as a, 
b, and c, respectively. By construction, a, b, c can differ 
by at most I since vertices X, Y and Z are adjacent. 
Then the classification function class(v) of a node v is 
defined to be 



class(v) 



( Pi 


for a < b — c, 


Pi 


for b < a = c, 


p 3 


for c < a = b, 


< Pi 


for a = c < b, 


p 5 


for a — b < c, 


p e 


for b = c < a, 


{ P7 


for a = b = c. 



(4) 



It should be mentioned that the definition of node clas- 
sification is recursive. For instance, Class P\ and P4 in 
A^ belongs to Class Pi in A t +i, Class Pi and Pq in 



A^ belongs to Class P2 in A t +i, Class P3 in A^ 



be- 
longs to Class Pj in A t +i, and Class P5 and Pj in A^ 
belongs to Class P5 in A t +\. Since the three vertices 
X , Y and Z are symmetrical, in Apollonian network we 
have the following equivalent relations from the view- 
point of class cardinality: Classes Pi, Pi, and P3 are 
equivalent to one another, and it is the same with Classes 
P4, Pg, and Pq. We denote the number of nodes in net- 
work A t that belong to Class Pi as Nt ^ , the number of 
nodes in Class Pi as N tj p 2 , and so on. By symmetry, we 
have N t , Pl = N,p 2 = and N t , Pi = N t , Pb = N t , Pe . 

Therefore in the following computation we will only con- 
sider N t .p l: A tj p 4 , and N t ,p 7 . It is easy to conclude that 

N t = N t , Pl + N t .p 2 + N t ,p 3 + N t , Pi 

+N t , Ps + N t}P6 + N t ,p 7 + 3 
= 3Ai, Pl +3At :Pl +A t ,p 7 +3. (5) 

Considering the self-similar structure of Apollonian net- 
work, we can easily know that at time t + the quan- 
tities Nt+i, Pl , A t+ i j p 4 , and N t +i y p 7 evolve according to 
the following recursive equations 



N t+ltPl = 2N t , Pl +2% 4J 
N t +i,p 4 = Nt, Pi +N t ,P 7 , 
N t+hP7 = 3N t , Pl +l, 



(G) 



where we have used the equivalence relations: Nt, Pl = 
Nt, P2 = N t , Pa and Nt,P 4 = N t .p 5 = N t , Pe . 

For a node v in network -At+i, we are also interested 
in the smallest value of the shortest path length from v 
to any of the three peripheral vertices X, Y , and Z. We 
denote the shortest distance as /„, which can be defined 
to be f v — min(a,b,c). Let dt,p 1 denote the sum of f v 
of all nodes belonging to Class Pi in network A t . Analo- 
gously, we can also define the quantities dt : p 2 , d t .p a , ■ ■ ■ , 
d tt p 7 . Again by symmetry, we have dt.p ± — d t ,p 2 — d ty p 3 , 
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d t ,p 4 = d t ,p 5 = d t ,p 6 , and d tiPl , d t , Pi , d t , Pl can be written 
recursively as follows: 




2d t , Pl +2d t , Pi , 
dt,Pi + dt,p r , 
3(d tjPl +N t ,p x ) + l, 



(7) 



Next we begin to determine the value of the crossing 
path length <d]' 2 . Using the definition of /„, Q]' 2 can be 
rewritten as 



E 



(fu + fv+Suv), (8) 



where S uv can be or 1 depending on whether the two 
paths corresponding to f u and /„ may meet at the same 
vertex (X or O). If the two paths meet at the same 
vertex, 6 UV = 0, 8 UV = 1 otherwise. According to the 
above classification of nodes, O]' 2 can be expressed as 



7 7 



®*' 2 = EE 
i=l j=l 



Nt-x,Pj dt-i,p t 



(9) 



where P, and Pj are node classes of A\ 1 } 1 and A^z 1 , re- 
spectively. The quantity 5p i p j may equal to cither or 
1, which are related to two kinds of paths: paths from 
nodes in Pi to the three peripheral vertices of and 
paths from nodes in Pj to the three peripheral vertices 



of A 



(2) 
t-1- 



If these two sorts of paths may meet at the 



same peripheral vertex, 5p i p j — 0; If they do not meet, 
Sp i p j — 1. Combining previous equations and results, we 
get the final expression for Gy 2 , 



eP = 



2178 



exp(-iTrf) 2*(-62 - 16iV2) 



+lli2*(10i + 7V2)exp^^ 

-lli2^(-10i + 7\/2)exp(^^ 
+i2 1+t (31 i + 8\/2) exp(2i7rf) + exp(i?rf) 
(-484 + 11 • 3 2+t + 9 3+t + 22 1 ■ 9 1+t ) J (10) 



Inserting Eq. (TIT))) into Eq. $2$ and using the initial 
condition D\ = 6, we have 



D, 



330(-3 + i%/2)(-3i + V2) 
-1551*25+1 + 5i(-i) f 2 4 +t + 5ii*2 4 +^ 



6655 i+ 155(-i) t 23+l 



+31i(-l) t 2 2+t + 133i3 3+t 
+55it3 1+t + 55ii3 1+2t 



40i3^ 



(11) 



□ Analytical Results 
• Numerical Results 



10000 100000 1000000 



A', 



FIG. 2: Average path length dt versus network order Nt on 
a semilogarithmic scale. The solid line is a guide to the eye, 
which clearly shows that the APL scales logarithmically with 
the network size. 



Then the exactly analytic expression for average path 
length can be obtained as 



* = 1815(15 + 8 .3' + 9') P 5 + 5H)2 ^ 

+5i*2 4+ 5 - 155i(-i) t 2 i * i + 155^*2^ 
+31(-l) t 2 2+t + 133 • 3 3+t + 40 • 9 1+t 

+55£ • 3 1+f (l + 3 f )^ , (12) 
which can be simplified according to t as follows: 

4[6655+5-2 5+ ^+31-2 2+ ' + 133-3 3+t +40-9 1+ '+55-3 1+t (l+3 t )t] 



1815(15+8-3*+9 t ) 
-31-2 2+t + 133-3 3+t +40-9 1+ '+55-3 1+t (l+3 f )t] 



4[6655-155-2~ 



d t = I 



4[6655- 



1815(15+8-3*+9 t ) 
f + 133-3 3+t +40-9 1+t +55-3 1+t (l+3')t] 



1815(15+8-3* +9*) 
-31-2 2+t + 133-3 3+t +40-9 1+t +55-3 1+t (l+3 t )t] 



4[6655+155-2~ 



1815(15+8-3*+9 t ) 



(13) 

for t = 0,1, 2,3 (mod 4) are given consecutively. In the 
infinite network size (i.e., large t), d t ~ jjt ~ lniV t . 
Thus, the APL grows logarithmically with increasing size 
of the network. We have checked our analytic result 
against numerical calculations for different network size 
up to t = 13 which corresponds to AT 13 = 797164. In all 
the cases we obtain a complete agreement between our 
theoretical formula and the results of numerical investi- 
gation, see Fig. 

We finally make a comment on the result presented 
in [15j. Using a combination of numerical experiments 
and heuristic arguments, Ref. [l5| advances claims about 
the interesting quantity of average path length. The au- 
thors of [TU claimed that d t cx (lniV t ) 3 / 4 . They thus 
concluded that ANs belong to a new class of networks 
which interpolates between small (d t cx lnN t ) and ul- 
trasmall (dt oc lnlnJVt) networks [7[. While the com- 
binatorial arguments and simulation results are sound, 
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the behavior of APL as t grows is quite difficult to esti- 
mate numerically. Just as it is hard to verify with direct 
computation that the harmonic series diverges, the ex- 
perimental results in [l~5| underestimate the asymptotic 
behavior of APL. In contrast to previous claims, our ob- 
tained precise value for APL shows that the conclusion 
of [15j ] seems questionable. The self-similar structure of 
Apollonian network allows us to compute precisely the 
quantity APL, which is difficult to apprehend through 
simulation. 

In conclusion, in this report we have derived analyti- 
cally exact solution for the average path length of Apol- 
lonian networks which has been attracting much research 
interest. We found that in the infinite network size limit 
Apollonian networks are small-world, the APL scales log- 
arithmically with network size. However, some authors 
have been proved that conventional random scale-free 
networks with degree exponent 7 < 3 show ultra small- 



world property @, BQ- In future, it is interesting to do 
further studies to reveal this dissimilarity between Apol- 
lonian networks and those stochastic scale-free networks, 
which have different scaling for APL. Finally, we believe 
that the analytical calculation here can guide and shed 
light on related studies for deterministic network models. 
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